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Abstract 

We consider the contour representation of the infinite volume Ising model at any 
fixed inverse temperature f3 > f3*, the solution of Yle-e^o e_ ^' e ' = 1- Let fi be the 
infinite-volume "+" measure. Fix V C Z d , A > and a (large) N such that calling 
G]yy the set of contours of length at least N intersecting V, there are in average 
A contours in G^ry under /i. We show that the total variation distance between 
the law of (7 : 7 G Gjyy) under \i and a Poisson process is bounded by a constant 
depending on (3 and A times e~^~^ ^ N . The proof builds on the Chen-Stein method 
as presented by Arratia, Goldstein and Gordon. The control of the correlations 
is obtained through the loss-network space-time representation of contours due to 
Fernandez, Ferrari and Garcia. 

Key words: Peierls contours; Animal models; Loss networks; Large contours; Ising 

model; Poisson approximation; Chen-Stein method 

AMS Classification: Primary: 60K35 82B 82C 60F17 60F05 



1 Introduction 

The infinite volume Ising model is one of the most studied in Statistical 
Physics. One of the fundamental historical fact was the Peierls argument for 
the existence of more than one phase at low enough temperature — which we 
assume throughout this paper — in any dimension d > 2. This argument was 
set by Dobrushin as the existence of at least two extremal Gibbs states, v + 
and v~ that are obtained as infinite volume limits of Gibbs measures with 
+ resp. — boundary conditions. In the Physics literature, a configuration of 
spins which is typical with respect to v + , is frequently described as a set of 
islands of — within a sea of + ; the contours being the boundaries of the is- 
lands. Various questions can be asked about typical configuration of contours 
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at low temperature. A very simple one that does not seem to have attracted 
attention is: what is the number of occurrences of contours larger than a fixed 
(big) N, that intersect a fixed volume V = V(N) say centered at the origin. 
To describe the physical idea of rarity of the island of — , a Poisson type 
behavior naturally comes in mind. We prove indeed that the distribution of 
these occurrences, under the condition that the mean number of occurrences 
in the volume V(N) is a fixed value A, approaches sharply a Poisson process 
of mean A. The paper is organized as follows: in section 2, we give the defi- 
nitions, state the main result and the Arratia-Goldstein-Gordon formulation 
of the Chen-Stein Method used in the proof. In section 3 we recall the graph- 
ical construction of the loss-network space-time representation of contours of 
Fernandez-Ferrari-Garcia which we apply to obtain the relevant bounds. In 
section 4 we prove the main result. 



2 Definitions and the main result 

Peierls introduced a map between typical configurations of u + or v~ into an 
ensemble of objects — the contours — interacting only by perimeter-exclusion. 
See, for instance, Section 5B of Dobrushin [6], for a concise and rigorous ac- 
count of this mapping. Contours are hyper surfaces formed by a finite number 
of (d — l)-dimensional unit cubes — links for d = 2, plaquettes for higher 
dimensions — centered at points of Z d and perpendicular to the edges of the 
dual lattice Z d + (|, • • • , |). To formalize their definition, let us call two plaque- 
ttes adjacent if they share a (d — 2)-dimensional face. A set of plaquettes, 7, 
is connected if for any two plaquettes in 7 there exists a sequence of adjacent 
plaquettes in 7 joining them. The set 7 is closed if every (d — 2)-dimensional 
face is covered by an even number of plaquettes in 7. Contours are connected 
and closed sets of plaquettes. For example, in two dimensions contours are 
closed polygonals. Two contours 7 and 9 are said compatible if no plaquette 
of 7 is adjacent to a plaquette of 9. In this case we write 7 ~ 9. In two di- 
mensions, therefore, contours are compatible if and only if they do not share 
the endpoint of a link. In three dimensions two compatible contours can share 
vertices, but not sides of plaquettes. Ising spin configurations in a bounded re- 
gion with "+" (or "— ") boundary condition are in one-to-one correspondence 
with families of pairwise compatible contours. 

The set of possible contours contained in A C Z rf will be denoted G[A]. The set 
of configurations of compatible contours is X[A] = {rj e {0, 1} G I A ] ; 77(7) rj(9) = 
if 7 rfj 9}. We denote I7I the number of plaquettes of the contour 7 and define 
the finite- volume Gibbs measure on X[A}: for r\ e X[A], 
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We fix (3 > /3*, the solution of a (/3) = 1, where a (/3) := Ee-.BBO e ~ m • Then 
there exists a unique weak limit \x = lim A ^ Z d fj,^ [7,8]. Notice that (3* is strictly 
bigger than j3p, the Peierls inverse-temperature, which is the infimum of the 
(3 making the previous sum finite. Let p 7 = / [i(drj) 77(7) be the probability 
of the presence of contour 7 under //. Fix iV > 0, let V = V(N) C Z ' 
and A > be such that if we define G = G(iV, V, A) := {k G G[Z d ] : \k\ > 
N and kC\V 7^ 0} as the set of contours with length at least N intersecting V, 
then Z) 7 eGP7 = ^- Hence V is a set such that in average there are A contours 
of length at least N intersecting V. Define the processes X := (77(7) : 7 G 
G), Y := (F(7) : 7 G G) where 77 has distribution fi and ^(7) are iid 
with Poisson distribution of mean p 7 . We denote £(X) the law of a process 
X and ||£(X) — £(Y)|| TV the total variation distance between the laws of the 
processes X and Y [2]. 

Theorem 1 Fix (3 > (3*, then for any (3' G (/?*, (3), if M((3,(3') is given in 
(26) below we have: 

||£(X) -£(Y)||tv < M((3,(3')N d+1 Xe-^-^ N (2) 



The proof of this theorem is an application of the Chen-Stein method [10,5,2,3] 
as proposed in [2]. In order to describe it, let the distance between contours 
be (1(^,6) := mm{\x-y\ : x G 7, y G 9}. For D = D(N) > to be fixed later 
we define B 7 = B^iV,/)) = {9 G G : d( 7 , 9) < D}, a "D-neighborhood" of 
7. Let p 7 e = / ji{drj)rj{j)rj{9) and define 6, = bi(N,D) by 



7eG 6»eB 7 7<eG 7 ^6»eB 7 

6 3 :=E E | E ^(7)-P7kW) = ^GG\B 7 )}| 

7GG 



Theorem 2 of [2] shows that the total variation distance between the Poisson 
process Y and the process X is dominated by 2(2&i + 26 2 + b 3 ). To prove 
Theorem 1 we choose D equal to iV times an appropriate constant and then 
find bounds on 6j. The bound on b\ does not offer problems, but the bound 
on 62 requires the inequality p 7 # < p 7 e -13 ^ +p e proved in Lemma 4 

below. The Chen-Stein method has been applied mostly to examples where 
63 = [2], while in our case b 3 7^ 0. To dominate 63 we control the correlations 
between "far" contours in Lemma 7 below. The proof of those lemmas is based 
on the graphical representation of a loss network, a Markov process having 
/i as invariant measure [7,8]. There are many papers dealing with Poisson 
approximations; we limit ourselves to quote the books [1] and [3] — with a 
nice introduction explaining the Chen-Stein method — , the original paper of 
Chen [5] and a monograph of Stein [10]. There are a few results for random 
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fields. The number of "— " spins in the "+" measure of the Ising model at low 
temperature and/or high external "+" magnetic field converges to Poisson 
[4,9]. In this case the variables are positively associated and the Chen-Stein 
method works via cluster expansions giving explicit convergence rates [9]. In 
[8] the approximation is established for the contours of the Ising model at low 
temperature and without magnetic field. The available bounds for the total 
variation distance are of the order of the inverse of the volume of the observed 
set V [9,8]. 



3 Graphical construction 

In this section we quickly review the construction and some results of [7] 
relevant to our proof. To each contour 7 E G(Z d ) we associate an independent 
(of everything) marked Poisson process iV 7 on R with rate e - ^ 7 L We call 
7^(7) 6 R, 7 E G(Z d ), the ordered time-events of iV 7 with the convention 
that T (7) < < 7\(7). For each occurrence time 7* (7) of the process iV 7 
we choose an independent mark 5* (7) exponentially distributed with mean 1. 
At the Poisson time-event 7* (7) a contour 7 appears and it lasts 5* (7) time 
units. The resultant object is the random family C = {{(7, 7* (7), ^(7)) : 
i 6 Z} : 7 6 G}. A marked point (7, T k {^), S k {^)) E C is identified with 
7 x [2^(7), Tk(j) + Ski'j)], the cylinder with basis 7, birth-time Tk(j) and 
lifetime ^(7). The life of the cylinder is the time interval [7^(7), 7^(7) + 
Ski'j)]- For a generic cylinder C = (7, t, s), we use the notation Basis (C) = 7, 
Birth (C) = t, Death (C) = t + s, Life (C) = [t, t + s]. In the sequel P and E 
are the probability and expectation in the space where C is defined. 

For t 6 R we define 6(7, C) = Ecec l{Basis (C) = 7, Life (C) 9 *}■ The 
above process, called the free network, is a product of independent stationary 
birth-and-death processes on N G ' zd l with t E R whose generator is given by 

A°f(0= E e-^l[/(£ + 5 7 )-/(0] + E e(7)[/(e-5 7 )-/(0]- 
7 eG[z d ] 7 eG[z d ] 

The invariant (and reversible) measure for this process is the product measure 
fi° on N G[zd] with Poisson marginals: /i°{£(7) = k} = exp (V^l) (e-^^^/A;! 
In particular, for any t E R, E/(£t) = At /. 

We say that cylinders C and C" are incompatible and write C" 7^ C if and 
only if Basis (C) 7^ Basis (C) and Life (C) n Life (C) ^ 0. We say that two 
sets of cylinders A and A' are incompatible and write A A' if there is a 
cylinder in A incompatible with a cylinder in A'. Otherwise we use the sign 
~ for compatibility. For any cylinder C define the set of ancestors of C as 
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the set of cylinders in C born before C that are incompatible with C: Af = 
{C G C; C rJL C, Birth (C) < Birth (C)}. Recursively for n > 2, the nth 
generation of ancestors of C is A% = {C" : C" G Af for some C G A^ r _ 1 }. 
Let the c/an of C be the union of its ancestors: A c = [j n >i &-n - Under the 
condition (3 > (3* all cylinders in C have a finite clan with probability one. This 
property called no backwards oriented percolation is essential to show that the 
loss network 774 can be constructed in a stationary way for tel. Assume that 
there is no backwards oriented percolation. The construction is as follows. 
All cylinders in C are classified as kept and erased. Since all clans are finite, 
we can write C = U„> C„, where C n := {C G C : A% ^ 0, A^ +1 = 0}. 
Inductively we set K = C , D = 0, K n = {C G C n \ U™T 1 (D i U K;) : 
{C} ~ U^Ki} and D n = C n \ [K n U U^D; U K;)]. Let the set of kept 
cylinders be K = U n K n and the set of erased cylinders be D = U n D„. Clearly 
KUD = C. The event {C G K} is measurable with respect to the sigma 
field generated by A c . In words, it is sufficient to know the (finite) clan of 
C to know if C is kept or erased. The stationary loss network is defined by 



774(7, C) = SceK l| Basis (C) = 7, Life (C) 9 t [. The process r) t is Markovian 



The unique invariant (and reversible) measure for this process is the Gibbs 
measure \i. In particular, for any t G R, Ef(r) t ) = /if. We study properties of 
by studying the law of 770, the stationary loss network at time zero. 

The presence/absence of contours intersecting a region A at time t depends 
only on the set A A,t , the union of the clans of the cylinders of C with ba- 
sis intersecting A and life containing time t, that is A A '* := {C G A c : 
Basis (C) PI A ^ 0, Life (C) 3 t}. In particular 774(7) is a (deterministic) func- 
tion of A 7 '* defined by 774(7,0) = 774(7, A 7,i ). The reason is that in order to 
determine if C G K it suffices to look at A c . When t = we will use the 
notation A A instead of A A '°. 

Let C ltt ,r = (7, —t,t+r) be a cylinder with Birth (C 7j t jT .) = —t, Basis (C ltt ,r) = 
7, Life (C ltt ,r) — t + r. Conditioning to the birth time of the (kept) cylinder 
alive at time zero with basis 7, 



with generator 



Af(v) = E Hv + £7 g x} [f( v + <g - /(T7)] 



7GG 



+ E*K7)[/fa-«g-/(0]- 



(3) 



7GG 



P 7 = 




•^ItI dtdr 



(4) 



(R+) 2 
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Coupling of clans. For disjoint sets A and T of Z d it is possible to construct 
(A A , A T , A A , A T ), a coupling between four sets of cylinders satisfying (a) 
A A = A A and A T = A T , in distribution; (b) A A U A T = A AuT ; (c) A> and 
A T are independent; (d) If A A ~ A T , then the marginals coincide: A A = 
A A and A T = A T . We write P and E for the probability and the expectation 
of the coupling. The following bound for the probability of incompatibility 
between clans follows as in the proof of (2.13) of [8]: for any f3' G ((3*, (3), 

P(A A rjL A T ) < 2(1- a ((3')y 2 E \ x ~ V\ e-V~^ x - y \ . (5) 



4 Proof of the theorem 

The proof of the theorem is based on a sequence of lemmata. 
Lemma 2 Denote p = a ((3). Then for 7 such that \ j\ > N, 

exp{-(/5 + p)| 7 |} < Pl < exp{-/3| 7 |} (6) 

PROOF. Bounding the conditional probability inside the integral in (4) by 
one we get the upper bound. For the lower bound notice that the event {C 1;t>r G 
K} is measurable with respect to the sigma field generated by the cylinders 
born before —t, the birth time of C 7jtjr . Hence 

P(C 7>t>r G K I C 7>t>r G C) = P( H {V-t(0, A e -*) = 0}) (7) 
>P( {e-t(e, C) = 0}) = exp(- e ~ m ) > e ~ pbl ( 8 ) 

The first inequality follows from ^(7) < £t(7); the second identity from the 
distribution of ^ described in Section 3. Inserting this inequality in (4) we get 
the left inequality in (6). 

Lemma 3 The following inequalities hold 

\V\ < A(]T ItI" 1 ^)" 1 < XNeW+W (9) 

7£Go 

PROOF. The first inequality is immediate. Bounding below the sum in (9) 
by one term (choose any 70 such that |7o| = N) we dominate the middle term 
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in (9) by \lo\/p 10 which using the left inequality in (6) is bounded by the rhs 
of (9). ' 



Lemma 4 We have p 7 <j < p 7 e ^ +p e e ^ 

PROOF. Notice that p^ = if 7 9. Recall the notation used in (4). 
Consider compatible 7 and 9 and condition to the birth times of the kept 
cylinders alive at time zero with bases 7 and 9 to obtain as in (4), 

Pje = J P(C 7; t ir - G K, Cg :S:W G K I Cjj,r G C, Ce, s , w G C) 

(R+) 4 

x e - ' e _s e _r e~ w e'^ e~ m dsdtdrdw (10) 

For s < t the event {C 7itjr G K} is measurable with respect to the sigma field 
generated by {C : Birth (C) < —t}. Hence it is independent of {Cg jS:W G C} 
and we have 



l{s < t} P(C 7)t , r G K, C e , s , w G K | C 7)t>r G C, C e , s , w G C) 

< l{s < *} P(C 7)t , r G K I C 7it , r G C, C e , s , w G C) 
= l{s < *} P(C 7it>r G K I C 7it , r G C) 

< P(C 7)t , r G K I C 7Ar G C) ' (11) 

Analogously l{s > t} F(C 1;t>r G K, C 0t9>w G K | C 7)t , r G C, G C) < 

P(Ce,s,w G K I Ce )S)U , G C). This and (11) allow us to factorize the integrals in 
(10) and then use (4) to get the lemma. 



Lemma 5 The following bound holds for b\, for all e < d{(5 + p) 



6i<EP + l7l) d P 7 E (12) 
7eG oeG 



PROOF. Inequality (12) follows immediately from (6). To prove (13) it is 
sufficient to show that for all e < d(/3 + p), 

+ < 2\N(j±^)(^±±) d eMeN). (14) 



7eG 
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Using (D + |7|) d < (D + l) d e^ £d /e d and (9), the lhs of (14) is bounded by 

V e l7|ed p 

Me.^,, ^ 1 ''^. 1 (15) 

tgGq 



Divide the sum in the numerator in two pieces. The first one contains those 7 
such that N < \^\ < aN with a — ((3 + p)/((3 — (3'). This part is bounded by 
aNXexp(deaN). It is not difficult to see that if ed < (3 — (3', the second part 
is also bounded by aNX exp(deaN). 



Lemma 6 The following bound holds for b 2 

h < A {B±l)' [2N (j±L) + .-<"->» (16) 



PROOF. By Lemma 4, 6 2 is bounded above by 

J2(D + \l\) d P, E + EP+l7l) d e-^l E» (17) 

7 eG eeG 7 eG 6»eG 

The first term was controlled in the previous lemma. The second one can be 
treated in a similar way. 



Lemma 7 The following bound holds for b 3 . For any f3' G ((3*, (3), letting 
f/ = ao(P% 

h < 2A e~^ N + Q(f3, (3') N XD d e A ^ N e~^ D (18) 

where Q{(3, (3') =4(1- p')- 2 {(3 - (3')~ 2 {2d - l) 2 {(3 + p) and A{(3, (3') = 
((/? + p)(P- [3T 1 + W + P) log(2rf - 1). 

PROOF. Since <7 (770(0, C) : 6 G G \ B 7 ) C a{A 9 ' : 6 G G \ B 7 )) 



63 < Y, E \ E {vo(l,A^)-p,\a(A e : 0gG\B 7 )} 
7eG 

= J2 E \ E {vo(l,A^) - m (^A")\ a(A e : 0gG\B 7 )}| (19) 



7€G 



<E#V U A*} (770(7, A 7 ) + 770(7, A 7 ))) (20) 

7€G 6>€G\B 7 
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where the identity comes from property (c) of the coupling defined in Section 
2 and the second inequality from property (d). We divide the sum in (20) in 
two parts. The first one is over those 7 G G of length bigger than KN, for 
some K to be chosen later. Dominating the indicator function by one, this 
sum is bounded by 

]T 2p 7 < 2A e^ N e-W-W (21) 

7 eG,| 7 |>ifAr 

The second sum is over those 7 G G with length less than KN. Dominating 
770 (•) by one, the summand is bounded by four times 



P(AV (J A e )<(l-p'r 2 J2 E \x-y\e-W>\ x -y\ 

0eG\B 7 (AT) x£ 7 y.d(j,{y})>D 

using (5) in the first inequality. Hence the second sum is bounded by 



4A {KN) {2d-l) NK+2 D d 
(1 - p'f W ~ W) 



(P-P')D Jf3+p)N 



e ^ 1 ' e 



(23) 



using (9) to bound \V\. The sum (20) is then bounded above by the sum of 
(21) and (23). Fixing K = |±£ + 1, we get (18). 



Proof of Theorem 1. It remains to choose D in such a way that 2(2b\ + 
2&2 + ^3) is bounded above by the right hand side of (2). We already know 
that bi and 62 decay as e~( /3_/3 * _e ^. If we choose D = SN in (18), where 



|±£[l + Iog( M -l)(^ + ^)] + l, (24) 



then 



b 3 < ^ + P) dNd+1 e -(^')" (25) 



With this bound and the choice e = f3' — (3* in (13) and (16) we get (2). The 
constant M(f3 ) f3') is given by 

M(R ft\ 1 n (g ± et ( dd + 2 aog(2^-l)) 2 \ 
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